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Vector-Like (VL) quarks arise in the main alternatives to the supersymmetric extensions of the 
Standard Model (SM). Given the experimental possibility of a 125 GeV Higgs boson with rates 
significantly different from the SM expectations, it is motivating to study the effects of VL quarks 
on the Higgs boson cross sections and branching ratios. We perforin a systematic search for the 
minimal field contents and gauge group representations of VL quarks able to significantly improve 
the fit of the measured Higgs rates, and simultaneously, to satisfy the direct constraints on VL 
quark masses as well as the electro-weak precision tests. In particular, large enhancements can be 
achieved in certain diphoton channels - as pointed out by both the ATLAS and CMS Collaborations 
- optimizing then the Higgs rate fit. This is a consequence of the introduction of VL quarks, with high 
electric charges of 8/3 or —7/3, which are exchanged in the Higgs-to-diphoton loop. Interestingly, 
the field contents and formal Higgs couplings obtained here are similar to those of scenarios in 
warped/composite frameworks arising from different motivations. The various exotic-charge quarks 
predicted, possibly below the TeV scale, might lead to a rich phenomenology soon at the LHC. 

I. INTRODUCTION 

The main drawback of the Standard Model (SM) is probably the gauge hierarchy problem induced by the divergent 
quantum corrections to the mass of the Higgs boson [IHI]. Various alternatives to the supersymmetry, from the 
AdS/CFT correspondence [5] paradigm, address this problem: the so-called little Higgs models [SHE], the composite 
Higgs [9Tfl"3"] and composite top [H] scenarios, the Gauge-Higgs unification mechanism (see e.g. Ref. [THUS]) or the 
warped extra-dimension setup proposed by L.Randall and R.Sundrum (RS) [TTl HB] _ and its well-motivated version 
with SM fields in the bulk [19] allowing to generate the fermion mass hierarchy [2"0Tl36| . 

All these frameworks predict the existence of additional color-triplet states with vector-like gauge couplings; these 
Vector-Like (VL) quarks arise [100] as Kaluza-Klein (KK) excitations in the higher-dimensional scenarios, as excited 
resonances of the bounded states constituting the SM particles in composite models and essentially as partners of 
the top quark being promoted to a larger multiplet in the little Higgs context. The VL quarks are also prevalent 
as KK excitations of bulk fields in other higher-dimensional frameworks (as e.g. in Ref. |37j ) or as new components 
embedded into the simplest SU(5) representations in gauge coupling unification theories 38J. Hence VL quarks are 
predicted in most of the alternatives to supersymmetry and their masses can reach values below the TeV scale as it 
could occur e.g. for the so-called custodians, KK states arising in the custodially protected warped extra-dimension 
models [HH] [TUT] . 

Now from the experimental side, all the recent data on Higgs boson searches performed at the Tevatron and 
Large Hadron Collider (LHC) constitute an important insight in the exploration of the Electro- Weak Symmetry 
Breaking (EWSB) sector. Hints in the data collected in 2011 by the ATLAS and CMS experiments, which were 
observed in various final states, could correspond to slight excesses of events induced by a Higgs boson with a mass 
of nih — 125 GeV as announced in Ref. [45, . New and updated LHC results have appeared at the Moriond 2012 
conference; then both LHC experiments reported local ~ 3a excesses [151 HZ] but more analyzed event statistics is 
still needed to confirm that these excesses with respect to the background estimation are not statistical fluctuations. 
The CDF and DO Collaborations have also detected an excess around ~ 125 GeV in the bottom quark search channel, 
h —> bb, corresponding to a local significance of 2.7cr [48] . 

Under this hypothesis of the existence of a 125 GeV Higgs boson, deviations with respect to the SM are observed: both 
ATLAS and CMS report an increase of all the Higgs production rates in the diphoton channel (h — > 77) compared to 
the SM Higgs rate, being close to the 2a level in some cases [see the precise data in the following detailed discussion]. 
The combined CDF and DO analyses point towards a best-fit cross section, for pp — > hV — > bbV [V = Z°, W bosons], 
centered at ~ 2 times the SM cross section [15] also for = 125 GeV [102] . Finally, the Higgs production rate in 
the channel h -> W+W~ is below the SM rate by more than la (still at 125 GeV) in the CDF+D0, ATLAS and 
CMS results. The other measurements are in a reasonable agreement with the SM expectations. 

The investigation of the Higgs sector, in addition to shed light on the cornerstone of the SM, represents a window on 
physics beyond the SM. The Higgs properties may e.g. help to discriminate between various new theories that could 
manifest themselves indirectly at the TeV scale. In this spirit, it is interesting to try to explain the above deviations 
of the Higgs rates from their SM predictions on the basis of corrections to the Higgs couplings. This approach is 
particularly relevant in the contexts of the well motivated VL quarks, discussed above, which are able to modify the 
Higgs interactions through fermion mixing and new loop contributions. 

In this paper, we adopt the generic approach of considering VL quarks as possible manifestations of several classes 
of theories beyond the SM, without specifying those theories. We elaborate the three types of minimal VL quark 
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models (and optimize their parameters) allowing to correct the 125 GeV Higgs couplings so that (i) most of the Higgs 
rates are significantly closer to their central measured values than the SM expectations and in particular for the 
channels h —> 77,66, W + W~ , (ii) all the Higgs rates belong to the la experimental regions and in turn improve the 
global fit of the SM, (Hi) the direct limits on the masses of involved VL quarks are respected and their contributions 
to the EW precision observables are acceptable. 

If the presently observed deviations from the SM Higgs rates are indeed to be explained by some VL quarks then 
our results predict the existence of at least two b' states (= VL quarks with the same electric charge Q e . m . as the 
6) plus at least either a pair of q s / 3 (VL quark with Q e . m . = +8/3) accompanied by a q 5 / 3 or a pair of g-7/ 3 
with a 5-4/3 . t' components (same Q e . m . as the top) may arise too. Some realistic mass intervals obtained are [103 : 
m i w 650-1100 GeV, m a i ~ 800 — 1000 GeV, m a i « 650—1000 GeV, m i « 800 — 900 GeV, m b w 840 GeV 

?S/3 q 5/3 9-7/3 9-4/3 ' 2 

and m t2 ~ 900 — 1010 GeV [104] . We do not claim those values to be the only allowed ranges as such a conclusion 
would rely on exhaustive and thus long explorations of the parameter spaces which are quite large in the considered 
scenarios. Let us mention that for certain gauge group representations, the gs/3 or 9-7/3 quark can be stable. 
Interestingly, one should also remark that the essential field content and couplings addressing the Higgs rate 'anoma- 
lies', in all the models constructed here, also hold in specific set-ups of the custodially protected RS scenario (or 
composite model). In that sense, the present Higgs data could be seen as an indication for the RS scenario and 
these specific custodial set-ups. Note that this kind of custodial set-up [containing principally two 6' and two (?g/3 or 
q_7/3 custodians] has not been considered so far in the literature on RS (composite) model phenomenology [39j|4Q] - 
including the implications for Higgs searches at colliders [ill |4"9"H5"3"] with LHC inputs [51] - except in Ref. [53] ([ [55] 
for a dual composite analysis) - disconnected from Higgs searches - where it arises from some considerations on the 
EW Precision Tests (EWPT). 

The VL field configurations that we point out could be realized in the context of any theory where the effects on 
the Higgs sector with a different origin from the VL quarks are negligible, like e.g. in RS scenarios with custodians 
around the TeV scale and decoupling KK excitations of the gauge bosons much above ~ 3 TeV (which is generally the 
order of the bound from EWPT |39l 1561 157]): then the minimal realistic VL quark configurations and the optimized 
Higgs rates would be mainly identical as those obtained in the present simplified framework. The other possibility is 
to have complete theories underlying the SM where VL quarks appear, together with e.g. light extra gauge bosons 
or an extended Higgs sector, so that several types of effects could affect the SM Higgs rates, in which case our results 
could be useful for theorists as a guide on VL quark influences. 

At this stage, one should mention related works. In Ref. [55] (see also references therein), it is shown that specific 
scenarios addressing the naturalness problem - via the presence of top-partners such as t' - like the little Higgs, multi- 
Higgs and pseudo-Goldstone Higgs boson (possible pattern of composite Higgs models) scenarios, are constrained 
in a non-trivial way by recent Higgs rate estimations [even before the Tevatron data]. Constraints from the recent 
Higgs data have also been imposed explicitly on the Minimal Composite Higgs Models with fermions embedded in 
spinorial or fundamental representations of SO(5) 59 61 . Effective approaches constraining generic Higgs operators, 
including possible deviations to the EW gauge boson couplings (but no effects from additional fermions not mixed 
with the SM ones), can be found in Ref. [62] (see Ref. [63] for effective interactions involving both the top quark 
and Higgs field). Another study at the effective coupling level shows that ad- hoc massive fermions in color SU(3) C 
representations up to the 27 and with |Q e .m.| < 2 can explain the simultaneously observed inequalities (in the author 
notations): gV < 1, Vj > 2, gj > 1 [64 (66 data from Tevatron not included) 105 . Some of the recent LHC 
investigations on the Higgs scalar have also been considered in completely different theoretical contexts such as the 
Universal Extra Dimensions, where the constraints from the 77 and W + W~ channels leave only an allowed narrow 
window near nth = 125 GeV [65] . Within the minimal supersymmetric extensions of the SM, the corrections to the 
hW + W~ vertex due to the extended Higgs sector constitute also the main source of corrections to the loop-induced 
hyy coupling, a problematic correlation preventing to obtain opposite corrections to the W + W~ and 77 rates as 
wanted today [BB] (and references therein) . Finally, the fourth generation [57] , the radion [BSJ or the dilaton [M] have 
difficulties to interpret the rate enhancement observed recently at LHC, for m/, = 125 GeV, in the Vector Boson Fusion 
(VBF) channel: VV — > h — > 77, while a fcrmiophobic Higgs boson could explain it [69]. Increasing considerably this 
channel rate is also doable by introducing new states with large electric charges that can be exchanged in the /177-loop, 
in the spirit of the highly- charged quarks introduced here; this idea was realized in the different context of the type 
II see-saw mechanism where doubly-charged Higgs scalars arise [70] [106] . 

In Section [H] we discuss the effects of different models of VL quarks on the Higgs boson observables. Then these 
models and their parameter space are confronted to the collider data in Section |III| We conclude in Section |IV] 

II. MODEL BUILDING 
A. The Higgs boson data 

We first define all the Higgs rate observables which have been measured at the Tevatron and LHC assuming a 
125 GeV Higgs boson. Those are the following signal strength modifiers (given with the last references for their 
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experimental value): fx hl = a h B h ^ / Vf M B™^ [H [HJ I7J, fty = oiflh^wK M BKw Ei EU [73], = 
a h B h ^ TT /al M Bl™ TT [Hi EH], = a hV B h ^ bh /a™B™ bb gSJ ED EU |7J , n qW = ^ q B^ W w/< q Cww [22] 
and n ql — o'hqq-Bh-i.77/o'™i?™ 77 [HI [75], where B stands for the branching ratios, ah for the total cross section of 
the Higgs production (dominated by the gluon-gluon fusion mechanism gg — > h), a^v for the cross section of the Higgs 
production in association with a V-boson and Oh qq for the VBF rate. One needs to introduce also the quantity [551 176] . 

0"hX -Bh_>. 77 _ 0.3 X CTgg^h + CThqq + 0"hZ + ^hW #h-»77 



7 hX ^11^77 U - 6 X °gg-+h + °hqq + CT hZ + °hW -°h-)-77 



where the factor 0.3 has been estimated recently from simulating additional QCD jets [55] to account for the efficiency 
of events issued from the gluon-gluon fusion to pass the cuts for the selection of the hV and hqq productions (basically 
on the Higgs boson transverse momentum). Similarly, to be rigorous there is a factor of about 0.033 for the suppression 
of the gluon-gluon fusion events (containing jets at NLO) by the dijet-class tagging in the /i g7 measurement; the 
uncertainty on this factor is of 70% [77] but it does not alter significantly anyway the theoretical estimation of fi ql in 
our framework (given the cr gg _>.h corrections and the absence of <J} lqq modifications). 

The experimental values for all these /i's observables are synthesized in the right-part of Fig.Q (or equivalently 
of Fig.Q). In order to summarize the various results on these figures, we have combined under the basic gaussian 
assumption the ATLAS and CMS data for every search channel where both experiments provide results - except for the 
h — > bb channel where wc find it more instructive to discuss the data separately (see next paragraph). The combination 
has been done without including the correlation; this is correct for the statistical error and the uncorrelated systematic 
errors while the correlated systematic ones, like the theoretical uncertainty, are expected to be subleading compared 
to all the others [55] (so the way those are combined should not be crucial). 

What comes out at a first glance on the experimental results of Fig. ([I]) is, in particular, the enhancement of the 
estimated rates for the three diphoton channels compared to their respective SM predictions. One notices also the 
significant reduction of dhBh-jww relatively to its SM expectation, observed simultaneously by the Tevatron and 
LHC. Concerning the h — » bb channel, if one does not consider the ATLAS best-fit value which is negative, an 
enhancement appears with respect to the SM especially at Tevatron where the obtained accuracy is better. The other 
channels fall into the lcr regions. 



B. The VL quark effects 

In order to improve the fit to the Higgs data, one would first need to increase the bb channel in Fig.([l]) [right-part] 
with respect to the SM. For that purpose, one needs to enhance the h — > bb branching fraction since the presence of VL 
quarks does not induce large tree-level corrections to the hVV vertex nor to the initial Vqq coupling involved in the 
hV production. To increase the Lh->bb width at rrih = 125 GeV, at least two VL &'-like states (say b' and b"), mixing 
together and with the SM b state, need to be introduced so that the absolute value of the bottom Yukawa coupling can 
be increased thanks to additional elements in the Yukawa coupling matrix arising in the (6, b' , b") basis |107| . At least 
a Yukawa coupling mixing b and b' and another one inducing the 6' — b" mixing are necessary; since the Higgs field is 
in a SU(2)l doublet, gauge invariance imposes that b' and b" must be embedded in different SU(2)l representations, 
as well as for the Left-handed 6l state (or the bn) and b'. Restricting to two 6'-like states and to multiplets smaller 
than a triplet to have the minimal field content set-up (extended set-ups do not bring different kinds of effects), one 
must thus embed b' in a singlet and b" in a doublet. 

Similarly, to increase theoretically the diphoton channels in the hqq and hX productions - as suggested by Fig.Q 
- we have to enhance Bh->77- One cannot increase f7 gg _j.h, to favor o^x, since CTh-Bh-s-77 is already significantly 
increased (from the i?h-j-7 7 enhancement required by the Ohqq-Bh-^ data) so that in view of the data cr gg ->.h must 
decrease for compensating this large crh-Bh-s.77 variation - as it will occur in this paper. Because of this aspect, the 
large experimental values of (ThqqSh-i.77 as well as a^xBh^-yy and the above Th-^bb behavior, the needed theoretical 
enhancement of i?h->77 via Th-*^ has to be important. 

The possible increase of the loop-induced /177 coupling [78] from the top quark loop-contribution, through a t — t' 
mixing increasing largely the top Yukawa coupling, is not possible as this SM top coupling is already close to its 
perturbativity bound and moreover the dominant (and opposite-sign) triangular loop-contribution to /177 remains 
to be from the W -boson exchange. Suppressing the top Yukawa coupling is not a solution neither, for increasing 
largely the ft.77 coupling. 

The last promising way to increase the /177 coupling is to introduce VL quarks with high electric charges leading 
to new loop-contributions favored by the two photon couplings. The two first possible exotic charges obtainable by 
extending SU(2)l multiplets of b' 7 t' are Q e . m . = —4/3 and Q e . m . — 5/3, as induced by the relation Y = Q e . m . — I3L 
(Y = hypercharge, I 3 l = SU(2) L isospin), and have in turn absolute values higher than the top quark one. However, 
given the present direct bounds around 600 GeV [see below] for the masses of such VL quarks 9^4/3, 95/3 (basically 
decaying like t' , b' respectively: <z_4/3 — > bW, 95/3 — > tW), we have found that their loop-contributions to the /177 
coupling do not reach large enough amounts in regard to Higgs fit improvements. Then one has to introduce the next- 
higher absolute charges, held by (7-7/3 and g§/3, to increase again the electromagnetic couplings of the ft.77 loop. The 
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Yukawa couplings and masses of these (7-7/3, (fe/3 must be such that their loop-amplitude interferes constructively 
with the W -boson exchange to generate enhancements. The condition for, say, the q a / 3 to be exchanged in the 
/177 loop is clearly that it must couple directly to the Higgs boson; this means that there should be at least two 
q s / 3 components, noted q s / 3 and g 8 / 3 , belonging to different gauge representations. To be minimal in terms of field 
content and without loss of generality, we restrict ourselves to two g 8 / 3 (or two (7-7/3) components embedded in 
SU(2)l representations up to triplets - including those reveals another type of model with respect to the 93/3/9-7/3 
decay. Besides, we do not consider charges, \Q e . m .\ > 8/3, as those do not bring effects of different nature and are 
less usual charges (even if Q e . m . = —10/3 and 11/3 are considered e.g. in Ref. [55]). Therefore, the only possibilities 
are to embed either g 8 / 3 in a singlet and q' s ^ 3 in a doublet or q s / 3 in a doublet and q' s ^ 3 in a triplet (or similarly for 

5-7/3)- 

A nice feature about the presence of highly-charged quarks is to greatly increase rh->. 77 and only this width, since 
the diphoton channels 'suffer' from some of the largest experimental discrepancies with the SM. 

As another good consequence of the field configurations selected above, ffh-Bh->ww wn l be significantly reduced as 
it seems indeed to be indicated by the data (see Fig.Q): i?h-s-ww is reduced due to the Th-j-bb increase and <7h due 
to the destructive interference between the g 8 /3 (or (7-7/3) loop and the top quark loop contributing to <7 gg ->.h [75]. Of 
course, this latter feature of destructive interference with the top contribution |108j must be preserved in the presence 
of additional VL t' or (7-4/3 for instance. 



The minimal models 



The combined theoretical conditions (discussed in Section II B ) for improving the fit of the Higgs data (presented 



in Section II A I lead to the following exhaustive list of minimal models for the VL quarks. A first class of models, 



denoted as Models of type I, is defined by the following four possibilities for the field content: 

(98/3,95/3)13/6 > (98/3)8/3 , (0-1/3 , (t',b")[ /6 or (6",9-4/ 3 )- 5 /6 > (!) 

(9-4/3,9-7/3)-ll/ 6 , (9-7/ 3 )-7/3 , (0-1/3 , (0 &")l/6 OT i b " > 9-4/ 3 )-5/6 . (2) 
(98/3,95/3 , 0*5/3 . (98/3,95/3)13/6 . (0-1/3 . (0 &")l/6 OT (6", <7-4/ 3 )l 5 /6 . ( 3 ) 
(0 9-4/3, 9-7/3) t -4/3 , (9-4/3, 9-7/3)-ll/6 > (&")-l/3 . (0 Ol/6 > (4) 

where we have written the field components in their transposed SU(2)l group representations together with the 
hypercharge as a global subscript. The Models I.Q (i.e. defined by the field content of Eq.Q) and I.Q are 
characterized by a stable 9i 7 / 3 in the case where gi 7 / 3 is the lightest field of all its multiplet partners (namely q 1 _ i / 3 
and 62). Indeed, the only potential decay channel in I.j2|, 9^7/3 — > q-4/ 3 W~, would then be kinematically closed. In 
I. pi, the decay channel through a virtual intermediate state, q 1 _ 7 / 3 — > Q^ai 3 W~ — > b^W~W~ , would be forbidden 



by the absence of b — b' mixing (recall that b = SM bottom quark). In I.Q and I.Q, open decay channels for the 
g_7/3 partners could be, gi 4 / 3 — > q 1 _ Y ^W + and b 2 — > q^J s W + — > q 1 _ 7 / 3 W + W + , where gi 7 / 3 would then appear as 
missing energy at colliders. 

The comparable case where gi 4 / 3 or 6 2 is the lightest field among all its multiplet partners (leading to a stable gi 4 / 3 
or 62) does not occur in the parameter space we will consider. A similar discussion hold for g 8 / 3 and its partners 
within the Models I.Q and I.Q. 

A second class of models, Models II, is defined by these two possible field contents: 

(98/3, 95/3, O5/3 ' (98/3,95/3)l3/6 ' (95/3, *")7/6 > (0-1/3 , (*'", Ol/6 0r (O 9-4/ 3 )- 5 /6 ■ ( 5 ) 

These models are characterized by the dominant decay channel of the highest-charge component: q^ 3 — > ql^* 3 W + — > 
tiW + W + (kinematically open in realistic frameworks), as allowed by the t — t' mixing for which (q'^ 3 ,t") has been 
added. The other possible decay into the ti instead of the t\ eigenstate is subleading due to the phase space suppression 
induced by the ti mass - or the three-body nature of the 95/3 decay if t 2 is virtual. 
The last type of models, Models III, is defined by, 

(6',9-4/3,9-7/3) t -4/3 , (9-4/3, 9-7/3)-n/6 , (&", 4U/z)U/* . ( & "')-i/3 and/or (t',b"% 6 . (6) 

Here the dominant decay of the highest-charge component is gi 7 / 3 — > 9^4/3 W~ — > b\W~W~ (kinematically open), 
as induced by the b—b' mixing - allowed by the presence of (b", g" 4 / 3 ) |109j . The similar decay obtained by replacing 
b\ with & 2 has a much smaller rate. 
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At this stage, it is interesting to realize a certain theoretical consistence: all the minimal models obtained here are 
similar to concrete warped extra-dimension [43] and their dual composite Higgs [55] scenarios (constructed to satisfy 
EWPT), in the sense that these concrete scenarios also possess the above crucial features allowing to improve the 
Higgs rate fit. Indeed, the representations I, II, III in Ref. [13] (with the extension of Eq.(16) therein) or B2 in Ref. [55] 
contain two (j_7/3 custodians coupled via a Yukawa term as well as two b' custodians mixed together and with the 
b quark through Higgs interactions, reflecting thus perfectly the VL quark configuration of the present Models I.(j2j, 
I.(|4j), III. Furthermore, the embeddings IV of Ref. [J3] or T3 of Ref. [55] have two q s / 3 custodians with a Yukawa 
coupling as well as two b' custodians with the required mixings, exactly as for the quark set-ups here in Models I.Q, 
I.|3j), II. The additional fields and mixings arising in these concrete realizations (like g_4/3, 95/3 states or heavy KK 
towers) are not expected to perturb drastically the potential Higgs rate ameliorations, and on the contrary, could 
even add more freedom. Besides, considering here all in all a unique set of VL fields - not a replica per generation 
- corresponds to the assumption in Ref. [431 155j where typically the custodians for the first two quark (and three 
lepton) SM generations decouple. 

We end up this subsection by writing explicitly the Lagrangian for one of these models. With the field content 
in Eq.(|5| for the Model II with say the (b", <7-4/3)!_ 5 / 6 doublet, all the possible mass terms and Yukawa couplings 
appearing in the generic Lagrangian, invariant under the SU(3) C x SU(2)l x U(l)y gauge symmetry, are, 





R/L \ / R/L 



r>. ( t ) Hb% + n' ( 6 ) Hb' R + Y b " ( q ^ J rfb c R + r_ 1/3 ( ^ J L/R Htb R/L + m b' L b R + m' b' L b' R 



-4/3 1 f) (/) + ™J9) m +^f!H ($») Z Z 1 +H.c. 



9-4/. J A «-■./> )u 5/3 V '" J A '" Jn J A \ f / \ f 

(7) 

where H represents the SM Higgs doublet, L/R the fermion chiralities, the Vs dimensionless Yukawa coupling 
constants and the m's various VL quark masses. Let us remark that the Yg/3, Y5/3 and 3^_i/3 terms could each be 
split into two terms with different chirality configurations and coupling constants. A field redefinition rotating b c R 
and b' R allows to eliminate the m term without loss of generality. The Yukawa couplings for the first two up-quark 
generations are not written in the Lagrangian ^ as their mixings with the top-partners t' , t" should be much smaller 
than the t-t' ,t" mixing. Indeed, new heavy i'-like states are closer in mass to the top quark and the top is in general 
more intimately connected to the ultraviolet physics, like in warped/composite frameworks. Since the CKM mixing 
angles [75] are typically small, the first two up-quark flavors should essentially decouple from the sector t,t',t". A 
similar discussion hold including the down-quark sector and the 6', b" components [110J. 

III. FITTING THE HIGGS BOSON RATES 

A. The theoretical parameter space 

We first consider the Model II which is quite attractive. In the left-part of Fig.Q, we present a domain of the 
parameter spac e where all the theoretical values of the Higgs rates belong to the experimental la regions [in the sense 
of Section II A| which are shown on the right-part of Fig.Q. The Model II considered in this figure contains the 



(b" . ^-4/3) doublet of Eq.S and its fundamental parameters appear in Eq.(j7j). For this parameter space exploration, 
we have typically let the relative /ix 7 rate lying within a still acceptable 1.4er region, to take into account an uncertainty 



in the QCD simulation of the efficiency for the gluon-gluon fusion contribution ( c.f. Section II A I . As the theoretical 
fi's quantities are normalized to the SM prediction, the QCD corrections arising in the VL quark contribution should 
essentially compensate the QCD corrections of the SM rate. 

Within the domain of parameter space presented in the left-part of Fig. 0, we observe that all the VL quark masses are 
well above their strongest direct experimental constraints which are at most, > 611 GeV (with the conservative 
assumption Bb 2 ^ tl w = 1) [80] . m t2 > 560 GeV (again with Bt^b^yv = 1) [SI]: TO ^ /3 > 611 GeV (assuming 

B q i 3 -^t ± W — 1> like a 62 state, which is a good approximation as the channel q^ 3 — > t^W is subleading) [50] and 
m i 1 4/3 ^ ^® GeV (assuming similarly B q i — 1, as a £2 state, in a good approximation) |81) . There are no 

existing searches so far for a particle with the uncommon main decay, q^ 3 — > t{W Jr W + (see discussion after 
Eq.([5])); anyway its mass values are quite high as illustrates Fig.([T]) - we have taken 650 GeV as the lower limit - and 
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FIG. 1: Left: Domain leading to 125 GeV Higgs boson rates inside the experimental la intervals for the Model II, with 
the (6",<j_4/3) doublet (c.f. Eq. |7j)), in the plan mg/3 versus m' s / 3 (in GeV). The values of the other parameters are fixed at 



Y = 1.01, Y 1 = 1, Yg/a = 2.5, Y 5/3 = -0.5, Y b = -0.053, Y b ' = 1, Y 6 " = 1, F_i /3 = 1, m' = 1200 GeV, m_ 4/3 = 900 GeV 



/3 
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7715/3 = 1000 GeV. Contour-level curves for the physical masses Tn q i / and are also shown. The other mass eigenvalues 

are almost constant over the shown plan like, mb 1 
900 GeV, or mainly depending on m.g/3 like, mt 2 
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-4/3 



900 — 1010 GeV in the domain shown, or around 2 TeV in this domain: 
/;/ =s 1250 - 3000 GeV, m„ 3 » 1500 - 3000 GeV, m„ 2 w 1600 - 3000 GeV. Right: Central values and la error bars for the 

J 9 5/3 q 8/3 

strength modifiers Hh-y, l^hv, V-Ht, fivt, ^ q w, ,"97 an d ^x-, (defined in Section II A I measured by the experiments indicated in 
front, for m/, = 125 GeV. The various strength modifiers are indicated by the associated cross sections and branching ratios. 
The plus symbols mean that the experimental results are combined. In each case the SM prediction corresponds to fj, = 1 
leading to the global Xsm value written in the figure. The small black circles correspond to the theoretical predictions of the 
strength modifiers for the point of parameter space also indicated as a circle on the left-side plot. This parameter set leads to 
the oblique S, T and Higgs fit x 2 values indicated near the circle on the top-left part of the figure. The little black squares are 



associated to a second parameter set where only one of the parameter values is changed: Yj 



8/3 



2.2. 



even higher for 5g/ 3 . The q^, 3 mass eigenstate decays either like the or as — > q^, s Z, q^ 3 h leading to a final 
state which has not been searched so far. Identical considerations hold for the <zf/ 3 , &3 and t% eigenstates. 
Concerning the couplings, all the absolute values of the fundamental Yukawa parameters entering Eq.([7]) have been 
taken larger than 0.5 not to introduce new unexplained hierarchies with respect to the top Yukawa coupling, Y , whose 
amount is close to unity as in the SM. The fundamental input parameter that is the bottom Yukawa coupling, Yf,, 
has also the same order of magnitude |lllj as in the SM (even if the absolute physical coupling is slightly enhanced 
typically by 6 — 6' mixings to increase the bb decay channel). The absolute Yukawa couplings in the mass basis do not 
exceed 2.5 and are thus below the usual perturbativity upper bound at v4tt- 

The acceptable domain in the left-part of Fig.|l]) is typically bounded from below by the la constraint on fi^z and 
from above by the condition on jix-y ■ This behavior of the Higgs boson rates is essentially due to the decoupling limit 
where mg/3 and Wg/ 3 tend to high values in which the rates tend to their SM predictions. Similarly, on the figure, 
the independence of the smallest mass eigenvalue m q i from m s / 3 , at high values of the latter relatively to m^ 3 , is 
explained by the decoupling effect of m s / 3 in the mass matrix [and reciprocally for m' 8 ^ 3 ]. It is also the case for m q i 
but not exactly for m q 2 ^ as this mass matrix also involves mg/3. 

It is remarkable that the domain in Fig.i|l]), leading to Higgs rates in a good agreement with the present data, is 
relatively large. Similar domains arise for different values of the parameters which have been fixed for drawing this 
figure. 

Let us finish this subsection by discussing the indirect constraints on the VL quarks. For the third generation quark 
sector, the tree-level corrections induced by the t — t' (b — b') mixings on the t (b) vertex are expected to dominate 
over the loop-level oblique corrections to the gauge boson propagators. Because of the relative heaviness of t' states, 
the predicted value for the Vtt CKM matrix element, including the t — if mixings, agrees with the experimental 
measurement obtained (without assuming 3x3 unitarity) through the single top production study [79] ■ In relation 
with the Zbb vertex, one could also try to address the LEP anomaly on the Forward-Backward asymmetry for the 
bottom quark as done in the specific RS context [4"Tl 02, 82 , 83J , assuming a discrepancy not due to under- estimated 
experimental errors, but this is beyond our scope. 

Concerning the interactions of leptons and first generations of quarks, one has to compute the corrections to the 
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gauge boson vacuum polarizations induced at one loop by exchanged VL quarks jHH 122] m the present model. The 
values of the oblique parameters S, T |112j that we can reach belong to the la regions induced by the long list of 
EW precision observables [113] measured mainly at LEP 79 j. This is true in particular for the parameters inside 
the domain of Fig. (JTJ) typically down to m 8 / 3 ~ 1900 GeV. Moving down to m 8 / 3 ~ 1000 GeV along this domain 
leads to T values up to ~ 0.5 (with S ~ 0.1 — 0.2) which would need to be compensated by other new physics effects 
than the VL quark ones. Such effects might be induced by new fields heavier than the VL quarks so that the former 
would correct T (which is extremely sensitive to new physics effects via EW observables measured typically at the 
per mille level) but would leave the quality of the improved fit to Higgs rates mainly unaffected (given its present 
large error bars typically of several ten's of percents). Such a scenario could be realized e.g. in a warped framework 
with relatively light custodians, heavy KK fermionic towers, heavy KK gauge bosons and possibly a gauge custodial 
symmetry (protecting the T parameter) in the bulk. 



B. The fits versus oblique parameters 

In the right-part of Fig.Q, we show the theoretical predictions of the Higgs boson rates still within the Model II 
for a point of parameter space first optimizing the fit - for a same fit quality then choosing the smallest S,T — (black 
circles) and for another point favoring the oblique parameters S, T (black squares) . Both points are compatible with 
the direct constraints on VL quark masses. 

The circle-points in this figure show that it is possible to obtain Higgs rates being all within the la regions. In 
this case, the fit is optimized in the following sense: starting from this situation where the four predictions for /i/ l7 , 
fJ-hz, fi>vb & n d /J-Xj are at the extreme la distances, one cannot improve one of these four quantities without moving 
another one of those out of its la region. One way of seeing this is as follows; the only possibility to increase ^x-y 
while keeping /i^ 7 at la is to increase i?h-s- 77 and decrease ov Now this <7h decrease would worsen the fit on fi^z 
as the only possible compensation by a i?h->zz increase via a Ph-s-bb decrease is forbidden if fiyt is to stay at the la 
level. 

The conclusion of this feature is that other parameters reaching the same quality of Higgs rate fit as above can be 
found but there exist no parameters improving the fit by comparison with the optimized situation described in the 
previous paragraph. Furthermore, neither different/additional SU(2)l multiplets nor higher electric charges of VL 
quarks - relatively to the present minimal model - could improve the fit in that sense |114j . These conclusions are the 
consequences of a certain tension among the Higgs rate data which restricts a little bit the potential fit ameliorations 
brought by VL quarks (whatever are the model and parameters). 

Let us also mention here that improving the Higgs fit through an enhancement of -Bh-s. 77 with a simultaneous 
suppression of Oh - as induced here by VL quark effects - seems to be favored by the generic analysis at the level of 
cross sections and branching fractions |66j . From a theoretical point of view, improving the fit on n qi by increasing 
the hVV coupling - from other origins of effects than the VL quarks - is at the price of extending the Higgs sector [55] . 
However, modifying this vertex, possibly through a custodial symmetry breaking [75], is probably the only way to 
reproduce the low /ihw experimental value without affecting too much the perfect agreement on fi^z of the SM [c.f. 
Fig.©]. 

The square-points in Fig. ([I]) correspond to S, T values clearly inside the ellipse associated to the la domain for 
combined EWPT 79 J. These points improve the fits of all the Higgs rates compared to the SM, especially fivb, 
/ix 7 , jtt g7 , Hhw and with the exception of fih-y (same typical deviation from data as in the SM), /x^z (larger but still 
acceptable deviation). This configuration leads to a clear improvement of the global x 2 function compared to the SM, 
as written in the figure [for /Xh 7 , ^hw only the most precise value, i.e. from ATLAS+CMS, is included in \ 2 whereas for 
fiVb only CDF+D0] j!15j . In addition to improvements similar to those, the circle-points all belong to the la regions 
of the Higgs rates with a remarkable increase for /i g7 ; nevertheless, the associated S, T values are slightly outside the 
la ellipse [75]. There might however be other kinds of effects from the physics underlying the SM responsible for such 
necessary small T compensations (S value being not problematic here) which could even let essentially unchanged 
the Higgs fit induced by the VL quark effects - as described at the end of Section |III A| Therefore, a noticeable but 
acceptable tension appears between optimizing the Higgs rate fit by introducing highly-charged VL quarks and still 
respecting the EWPT. 



C. The other models 



An example of parameter domain leading to Higgs rates within the la regions is also shown for the Model I with 
the (6", 9-4/3) doublet [defined by Eq.([2j)] in the left-part of Fig. (|2j). The fundamental parameters are noted in a 
similar way as for the Model II discusseain previous subsection. 

In this case, the constraints on m q i ^ deserve some more attention. Indeed, the <?_ 7 /3 particle is stable for the reasons 



the non-mixed and 



exposed in Section II C that apply here due to the kinematical feature, m q i < m_7/ 3 (m_ 7/3 being the mass of 

— 7/3 
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thus eigenstate 9-4/3), visible in Fig.((2|). Now such a stable particle has never been searched at 
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FIG. 2: Left: Same as in Fig. (|T|) but for the Model I with the (6", 9-4/3) doublet [see Eq.|2|] in the plan m_ 7 / 3 versus 



-7/3 
-4/3 



(in GeV). The fixed parameters read as, Y — 1, YI7/3 = 3, Yj = —0.053, Y b ' = 1, 



= 900 GeV. Three contour-level curves for m„ 



1/3 



-7/3 

= 1, m' = 1200 GeV, 



are shown. Other masses are almost constant over the plan, 



m b . 



4 GeV, m_, 2 ~ 840 GeV, mt, 3 ~ 1290 GeV, mt 1 ~ 173 GeV, and the remaining ones are, m 2 



1750 - 2750 GeV 



[in the presented domain], m q i = mm(m_4/3, 771-7/3), m q i = moi(m_4y 3 , m_ 7 / 3 ) [as there is no mixing term]. Right: 

Same as in Fig.|TJ) for the experimental data but with theoretical predictions from the Model I [with (6", 9-4/3)]. The black 
squares are associated to the same input parameter values as for the black circles except that YI7/3 — 2.2. 



colliders so far, but one can try to extrapolate constraints on it from investigations on other Heavy Stable Charged 
Particles (HSCP) even if this is a non-trivial task; constraints have been imposed on long-lived supersymmetric 
partners - namely a gluino g, a tau-slepton Ti, a top-squark t\ - and the most stringent ones have been derived 
recently at the LHC [57] (see Ref. [5SJ for analog studies at the Tevatron). After the hadronization stage, a long-lived 
gluino should form a 'R-gluonball' gg with a probability of typically ~ 30% (based on simulations) and a color-singlct 
'meson' gqq (q = u,d,s quarks) with Q e . m . — ±1 also for a significant fraction, whereas a long-lived stop might 
form a meson t\q (q — d,s) with a fraction typically about ~ 50% |89j . To get an idea of a possible mass bound 
on the stable 5i 7 / 3 which will also form a Q e . m . — ±1 R-hadron |116j . more precisely a baryon q\ 7 ^ 3 uu, one could 

apply the mass-dependent limits on the g (ti) production rates obtained in Ref. [87] to the q l _ 7 /^ pair production 
cross section calculated at NNLO for the 7 TeV LHC. We find, m q i_ 7 g > 850 GeV (800 GeV). Nevertheless, there 
are various limitations to the validity of this extrapolation: the fraction of q_ 7 , 3 uu should be simulated specifically 

(e.g. with the PYTHIA or HERWIG simulators) and is certainly different from the gqq (tiq) fraction, while the spin 
configuration and/or color- multiplet also differ from each other. Furthermore, the nuclear interactions experienced in 
matter by R-hadrons, suffering from large uncertainties, may lead to charge exchange; a recent work |90j modeling 
the HSCP nuclear interactions favors a scenario where the majority of R-hadrons made of a gluino or a squark would 
emerge neutral in the muon detectors. Assuming this result extends to all the R-hadrons made of the 9^7/3 {q^_7/ 3 <]q, 

an d using again the exclusion limits on the gluino (stop) production rates [87] obtained now under this charge 

suppression hypothesis, we find, m„i > 750 GeV (800 GeV). As for the limit on the fi production rate, it ends up 

— 7/3 

at a ~ 500 GeV mass [87) . Hence, taking the mean of the above four estimated limits, we obtain an indicative bound, 
m q i > 800 GeV, that is represented over the domain drawn in the left-part of Fig.|2j); one sees that a large part of 

the domain is passing this indicative test. 

There also exist studies on long-lived charged massive particles outside colliders (see e.g. the reviews in Ref. [89] and 
Ref. [91] ) leading in particular to limits on their abundance in ordinary matter; these particles can bind to a nucleus 
forming anomalously (super)heavy isotopes which have been searched, or even fall onto/through oceans and lakes to 
form heavy water molecules. For example, constraints on astrophysical fluxes of those particles have also been looked 
at |117| . However, no dedicated analysis (outside colliders) has really been performed to constrain significantly the 
mass of stable color-triplet fermions with high fractional electric charges. 

The eigenstate q-4/3 decays as, q-4,/3 — > q\ 7 / 3 W + or q-4/3 — > q^y 3 W + — > q\ 7 ^ 3 hW + ,q\ 7 ^ 3 Z°W + (open in most of 
the considered domain). In this case the stable 9i 7 / 3 could also be searched at colliders as missing energy associated 
with boson production, but there were no investigation so far on these kinds of (7-4/3 decays within the present 
framework [the 9-4/3 mass, m„ 7 / 3 , has typically high values anyway as shown in Fig.(|2|)]. Similarly, the q^_?/ 3 decays, 
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9^-7/3 ^ ^—7/3^' ^—7/3^^ ^nd. 

9-7/3 ~^ 9-4/3^ ~^ Q~7/3W + W , lead to specific signatures not yet analyzed. The 
constraint for the unmixed </_ 4 / 3 eigenstate mass, m_4/3 > 560 GeV (since B q i i/3 ^biW — 1) as possibly for a t^) [81j . 
is also clearly respected. Finally, the 6-sector is exactly as in the Model II of previous subsection so that the positive 
conclusions for the 62,3 mass bounds are similar. 

The square-points in the right-part of Fig. ^ for the considered Model I correspond to reasonable S, T values at 
the border of the la ellipse of EWPT [79]. Those theoretical predictions improve the fits of most of the Higgs rates 
compared to the SM, leading to a net improvement of the x 2 value. The circle-points even all belong to the la 
regions with again a remarkable increase for fi ql ; however, the associated S 1 , T values are now clearly outside the 
la ellipse [79]. Other types of new physics effects could reduce the T parameter down to acceptable values - the S 
parameter being already in a realistic range - as discussed in the end of Section |III A| 

The Models of type III [see Eq.([6])] lead to similar allowed domains of parameter space as above, as well as 
comparable qualities of the Higgs rate fits. The main difference with previous models concerns once more the type 
of bound applying to the mass of the VL quark with the highest absolute electric charge, namely here , 3 . As 

explained in Section II C it decays predominantly as, <?i 7 / 3 — > q^^W^ — > b{W~W~ , whose final state mimics that 
of, 62 - ► t{W~ — > biW + W~ . In obtained acceptable parameter ranges, the m q i values stand inside the interval 

— 7/3 

~ 700 — 1000 GeV being clearly above the experimental bound, m q i > 611 GeV (B q i ^^^ blWW ~ 1) [SO]. In the 

Models III, there are no stable VL quarks. 

Of course, we are not going to present detailed numerical results for the various minimal models of types I, II and 
III presented in Section |II C[ but similar qualitative and quantitative results as the ones presented throughout this 
paper hold. 



IV. CONCLUSIONS 



We have shown that the presence of VL quarks is sufficient to modify the SM Higgs boson rates such that all those 
belong to the present experimental la regions. The minimal field contents and gauge group representations of VL 
quarks allowing to achieve such improvements have been obtained. The key idea is to introduce VL quarks with 
high-enough electric charges to increase sufficiently the diphoton channel rates. Simultaneously, all the (in)direct 
constraints can be satisfied, even if we have pointed out a little tension between optimizing the Higgs rate fit and 
respecting the EWPT. The models obtained predict a rich phenomenology at LHC with several exotic-charge quarks 
possibly around the TeV. For example, the pair production of ?g/ 3 , decaying as q^ 3 — > tiW + W + , could lead to a 
spectacular final state with six W gauge bosons (two more W's than in the q^ 3 pair production [92]). 

One must mention that obviously the present measurements of the Higgs boson rates have large uncertainties. If 
the next data, in particular from the 8 TeV LHC, confirm the existence of a ~ 125 GeV Higgs field, two typical 
situations might arise. First, the main features of the present rate measurements, like diphoton channels significantly 
larger than in the SM, might be confirmed pointing towards manifestations of a theory beyond the SM. In this case, 
a simple adjustment of the present fundamental parameters should suffice to optimize the Higgs rate fit. Even the 
conclusion on the best-fit configuration underlined here would remain unchanged with respect to the central values; 
it is moreover true for any pure VL quark model (independently of its field content and gauge multiplets). A second 
possibility is that significant changes appear among the Higgs rate measurements so that the new diphoton rates get 
closer to their SM expectations. In such a case of agreement, our work could serve as a guide to constrain the VL 
quark models: for instance, such new data would disfavor the presence of 0(1) TeV VL quarks coupled to the Higgs 
boson and with absolute electric charges at 7/3 or above. 
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